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Abstract : We calculate the low-lying part of the spectrum of the Z 3 -symmetrical 
chiral Potts quantum chain in its self-dual and integrable versions, using numeri- 
cal diagonalisation of the hamiltonian for N < 12 sites and extrapolation iV — > oo. 
From the sequences of levels crossing we show that the massive phases have oscilla- 
tory correlation functions. We calculate the wave vector scaling exponent. In the 
high-temperature massive phase the pattern of the low-lying levels can be explained 
assuming the existence of two particles, with ^-charge Q = 1 and Q = 2, and their 
scattering states. In the superintegrable case the Q = 2-particle has twice the mass of 
the Q = 1-particle. Exponential convergence in N is observed for the single particle 
gaps, while power convergence is seen for the scattering levels. In the high tempera- 
ture limit of the self-dual model the parity violation in the particle dispersion relation 
is equivalent to the presence of a macroscopic momentum P m = ±y?/3, where (p is the 
chiral angle. 
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1 Introduction 



The ^-symmetrical chiral clock model ("CCVmodel") has been introduced in 1981 in 
order to describe the transition between commensurate and incommensurate phases 
of adsorbed gases on crystal surfaces jl], 0, |!|. It is a two-dimensional model defined 
by the partition function Z and the action S: 

Z = ^exp(-S) 

(a) 

S = -J2(g x cos[§vr(a E+li j / - a X}V )] + g y cos[|7r(a !ri2/+1 - a X}V ) - <pfj . (1.1) 

The Z3-spin variables at sites (x,y) take the values a x<y = 0, 1,2. g x and g y are the 
horizontal and vertical couplings, respectively. For non-zero chiral angle <p the system 
is parity-noninvariant in the y-direction. For p = 0, (|1.1|) is the usual Z 3 -clock model. 
A standard experimental example for a system exhibiting Z 3 -symmetry and which 
may be described by ( |1.1| ) is monolayer krypton covering of a graphite surface ||, |5| . 

Various techniques have been used in order to determine the phase diagram of 
the model: Mean field, Monte-Carlo, renormalization group, partial transfer matrix 
diagonalization, and hamiltonian limit finite-size calculations |7|-|11|. As yet, no full 



agreement has been reached between various authors on the details of the phase 
diagram. 

The CCVmodel is not self-dual. A self-dual model which should show commensu- 
rate-incommensurate phase transitions was proposed and studied by Howes, Kadanoff 
and Den Nijs in 1983 [[12]] . Their model is defined by the following quantum chain 
hamiltonian (or logarithm of the transfer matrix): 

2 N 

n = — -=Y, { e ~ ilp/3a i + e +ivV Vt + a (< ; " : 'T ; r;. l + e + ^ /3 r+r, +1 ) } , (1.2) 

V 3 j= i 

where Oj and Tj are 3 x 3-matrices acting at site j (the Z 3 -analogs of the Z 2 -Pauli 
matrices a z and a x ): 





T j= 10 (1.3) 



with to = exp(27rz/3). 7i contains the three parameters A, (f) and (p. A is the inverse 
temperature, and (p are two chiral angles. For <f> = (p the model is self-dual. 

For cf) = or (p = the hamiltonian TC follows from the CC3-action ( |1 . 1| ) ||7|. Q . 
but for general p and the corresponding two-dimensional lattice model has complex 
Boltzmann weights, so that there is no straightforward statistical mechanics interpre- 
tation. Because of some symmetries of TC, it suffices to consider only < 0, p < ir. 
The spectrum of Tt decomposes into three charge sectors Q = 0,1,2 because Ti com- 
mutes with the Z 3 charge operator Q = Ylf=i a j an d we ma y write the eigenvalues 
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of Q as exp(27u<5/3). Although for 0^0, parity is violated, TC is still translational 
invariant (we use periodic boundary conditions rv+i = I\) and the eigenstates will 
be labeled by the momentum eigenvalues p: 

p= -[f ],..., -1,0,1,..., [f]. (1.4) 

Using perturbation expansions and fermionic techniques, Howes et al. concluded 
that the phase diagram of the self-dual model (ip = 0) contains three phases: an 
ordered and a disordered phase separated by a Z 3 -Potts universality class line for 
low values of ip, and an incommensurate (IC)-phase for large ip. However, the most 
exciting result of their analysis (they used a lOth-order high-temperature expansion) 
was that for ip = 7r/2 the lowest gap of TC is a linear function of A (the 2nd to 10th- 
order coefficients in A were all found to vanish), just as it is the case for the Ising 
model. 

More precisely, let us define the gaps of TC with respect to the lowest Q = 0, p = Q- 
level, which we denote by E (Q = 0,p = 0). Observe that E (Q = 0,p = 0) is not always 
the ground state of TC, since the model has an IC-phase region where in general the 
ground state carries non-zero momentum p. 

By AEqj we denote the energy difference 

AE Q , = Ei{Q,p = 0) - E (Q = 0,p = 0) (1.5) 

where Ei(Q,p = 0) is the i th level (i = 0, 1, . . .) of the charge Q, zero-momentum 
sector. Then the Howes et al. result is 

AE 1>0 = 2(1 - A) forA < 1. (1.6) 



In 1984, v.Gehlen and Rittenberg |L3[] looked into the question, whether also other 
gaps of the Howes et al. model might be linear in A. They found that indeed this is 
the case. They also constructed a whole series of Z n -symmetrical models which have 
gaps linear in A. These models are defined by the quantum hamiltonians 

N n-1 

H = -J2J2{^^ + Xa k r} r?-*} (1.7) 

j=l k=l 

Here Oj and Tj are n x n-matrices acting at site j (the ^-generalizations of (|1.3|)) 
which satisfy 

Oj Y, y = I> Oj u) s s*\ a™ = T] = 1, uj = exp(2m/n). (1.8) 

The ^-version of ( |1.7| ) is the standard Ising quantum chain. 
The linearity of gaps is obtained if in (|1.7|) we choose 



ctk — cik — 1 — i cot(7rk/n). (1-9) 

At the same time eq. (|1.9|) guarantees |13[ that the models ( |1.7| ) are integrable in the 
sense that they satisfy the Dolan-Grady conditions fl4j| , which are equivalent |15| to 
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the Onsager algebra [T7|. 
to as " superintegrability 



m. 



This kind of integrability is now usually referred 
191 1201 ■ It ensures that if ak,ak satisfy eq.(|1.9|), then the 



A-dependence of all eigenvalues E(\) of H has the special Ising-like form Jl9|, [15 



E(X) = a + b\ + J2 4m iv / l + A 2 + 2Acosfl i 



1.10) 



Here a, b and 6j are real numbers and the mj take the values rrij 



-s, + l, 



with Sj a finite integer. For details on formula ( |1.10| ) see |T 

derived the hamiltonian f|l 



0- 



In 1989 Au-Yang et al. [21 



22j| derived the hamiltonian ( |1.7|) from a two-dimensio- 
nal lattice model built out of Z n -symmetrical (but not generally positive) Boltzmann 
weights. These satisfy a new kind of Yang-Baxter equations with the spectral param- 
eters defined on a Riemann surface of high genus. In a series of papers, the authors 
of [^TJ-p6| have been able, using in addition a functional equation |27fl , to calcu- 
late the complete spectrum of the Z 3 -hamiltonian ( |1.2[) in the superintegrable case 
= ip = tt/2. 

From the Yang-Baxter equations it follows that 
and dk satisfy the relations 



L7]) is also integrable if the a k 



a k = e l(2fe/n - 1) Vsin(7r£;/n), a k = e i(2k/n ~ 1)ip / sm(nk/n), (1.11) 



cosf 



A cos< 



TJT) for = ip 



;i.l2) 

: TT/2, 



Except for the superintegrable case, which is contained in 
these equations force H to be non-self-dual. Calculating the spectrum for the inte- 
grable, but not superintegrable case (|1.11|) , ( |1.12|) is a very difficult task. Nevertheless 
some initial steps have been done in attaccing this problem [P5I . 

In this talk, we shall give some new results on the Z 3 -model ( |1.2|) . We shall not 
rely on exact solutions and calculate low-lying energy levels of 7i by numerical diag- 
onalization (using the Lanczos algorithm) of the chain hamiltonian for up to iV = 12 
sites. The motivation for such a work is two-fold: 

1) For the superintegrable case, where exact results are available, the numerical evalu- 
ation of the exact formulae is not easy, and straightforward numerical diagonalization 
sometimes gives more transparent results on the behaviour of whole sets of levels. 
This may lead to the discovery of regularities which otherwise could be overlooked. 
Indeed, a good fraction of the progress in finding the special properties of ( p..7[ ) has 
been made possible by the extensive finite-size studies performed by the authors of 

B|. 

2) Of course, the main merit of the finite-size method is, that it is not restricted to the 
superintegrable case and so one can study how the special superintegrable spectrum 
emerges from the 92-dependence of the general spectrum. 

In the next Section we shall show that the massive phases are oscillatory, Sec. 3 
gives a calculation of the corresponding wave- vector scaling critical index, which sheds 

Ml 



new light on the scaling asymmetry noticed in [^, |24[. Sec. 4 then interprets the low- 
lying spectrum of the massive phases in terms of "elementary" particles and their 
composite states. Sec. 5 summarizes our results. 
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2 Oscillatory behaviour of the massive phases 



Due to its complex coefficients, the hamiltonian (|1.2|) allows ground state level cross- 
ings (the Perron- Frobenius theorem does not apply). The different phases can be 
easily distinguished by the ^-charge sectors and momenta of the lowest energy levels 
of 7i. In Fig. 1 for fixed number of sites iV = 10 we draw a map of the \/ip- plane (in 
order not to get lost on three-dimensional figures, in the following we shall consider 
the self-dual case <fi = tp if not explicitly stated otherwise). We label regions in this 
plane by the sectors of the two lowest levels, e.g. 0i/2 3 means that the ground state is 
in the Q = 0,p= 1-sector and the first excited state has Q = 2,p = 3. In most cases the 
lines mark alternately cross-overs in the ground state and in the first excited state. 



Fig. 1 : Map of the lines where cross-overs in the two lowest levels of the hamil- 
tonian fll.2|) (self-dual version ip = 0) for iV = 10 sites occur. The meaning of the 
labelling of the regions between the cross-overs is explained in the text. 
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Fig. 2: Schematic phase diagram of the self-dual version of the Z3-chiral Potts 
model defined by the hamiltonian eq.( |1.2| ). 

Inspection of Fig. 1 shows that there are at least 3 different phases (Fig. 2 gives a 
general schematic view of the phase diagram): 

• Low-temperature phase (I) at A > 1 and small 0: Here the three lowest states 
are from three different charge sectors and belong to the zero momentum sector. 
As A increases at fixed ip, the three levels form a kind of braid and alternate 
in forming the ground state. We mention that the gaps between the braiding 
levels decrease exponentially both with N and A, as it is expected for the low- 
temperature phase where the ground state is threefold degenerate in the thermo- 
dynamic limit N — > oo (the next levels follow at the distance of AE ~ 3(A — 1) 
as has been shown by Baxter |24|j ). 

In their low-92-branches the cross-over lines follow a simple pattern for A — > 00: 
all these lines tend to fixed values of (p according to the rule [|31|] (which we 
verified numerically): 
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Oo/lo 


(crossing in the ground state) 


> <p = 


71 /N 


Oo/2 


(crossing in the 1st excited state) = 


> <p = 


2n/N 


lo/2 


(ground state) 


> <p = 


3ir/N 


lo/Oo 


(1st excited state) 


> <p = 


Att/N 


2o/0 


(ground state) 


> <p = 


5tt/N 


2 /lo 


(1st excited state) 


> <p = 


67r/N etc. 



These rules are teaching us that for iV — > oo there will be an infinite supply of 
cross-over lines, which finally will fill the whole region A > 1 down to ip — +e . 
We have also verified that for free b.c. (and also for the integrable version of 
the model) the same asymptotic behaviour of the cross-over lines appears. So 
these cross-overs are not just artefacts of the boundary condition. We conclude 
that this phase is an oscillatory massive phase. In the next Section we shall give 
evidence that for increasing TV the cross-over lines converge to the line A = 1 
and give the calculation of the wave vector critical index. 

High-temperature phase (II). Dual to phase (I) this appears at A < 1 and small 
(f. Here the ground state always belongs to the sector Q = 0,p = so that there 
are no ground state level crossings. However, there is a regular pattern of cross- 
overs in the first excited state, which is always in Q — l, but with increasing N 
and (p switches to higher momenta. Perturbation theory in A teaches us that at 
A = the first excited state moves from momentum p to p\ at the chiral angle 

V = ^(Po + Pi). (2.1) 

Again, this rule ensures that for N —>■ oo there is an infinite supply of cross- 
over lines, this time in the first excited level. This gives rise to a modulation 
of Ap ^ 0-pair correlation functions (the Ap = 0-pair correlations are non- 
oscillatory) . 

Massless incommensurate phases (III, IV): In the middle upper part of the 
phase diagram there are plenty of ground- and first-excited state level crossings. 
Both the charge sectors and the momenta are changing simultaneously. With 
increasing N the pattern of cross-over lines becomes more dense. The loss of 
translational invariance of the ground state together with a A- dependent varia- 
tion of the ground state momentum is sign of the presence of an incommensurate 
(IC)-phase. 



Here the information from the exact solution of Albertini et al. 19 on the line 



ip = 7r/2 provides more detailed information than we are able to obtain from 



our chains of N < 12 sites: In JL{| it is shown that on the superintegrable line 
the ground state looses translational invariance in the interval 0.901293 < A < 
1/0.901293. In in addition it is shown that actually for A < 1 and A > 1 



there are two different IC-phases which do not join smoothly at A = 1. 
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The appearance of p ^0 in the first excited state in the high-temperature region 
is an indication that eventually the ground state may acquire momentum too, 
but iV < 12 sites are not enough to prove that this happens already on the line 
if = tt/2. 

Interesting information is obtained from our calculations on the different be- 
haviour of the self-dual (ip = 0) and integrable (cos <fi = A cos ip) hamiltonians 
for A > 1 and ip — > 7r: In the integrable case we find that the ground state 
cross-over lines for (p —>■ n concentrate below A ~ 1.6 and no cross-overs are 
left at (p — 7T above this value of A. This is in contrast to the self-dual case, 
where the cross-over lines for (p tc spread out to A = oo. In agreement with 
the phase diagram drawn by Roan and McCoy [p8]1 , we conclude that in the 
integrable case, phase (I) extends up to ip = ir and A ~ 1.6. 

For the self-dual case we are able to determine the boundary between phases (I) 
and (IV) with good precision, because the lowest gaps in region (I) are exponentially 
small as mentioned above, whereas the gaps in the IC-region are generally of order 
unity. In Tables 1 and 2 we give an idea of the precision with which the phase bound- 
ary to the IC-phases can be determined by our finite-size extrapolation [ pT|] . 









= 


ip (degrees 


) 






N 


A =0.1 


0.25 


0.5 


0.6 


0.7 


0.8 


0.9 


9 


170.860 


157.553 


135.605 


127.301 


119.986 


114.519 


112.333 


10 


171.112 


158.085 


135.767 


126.779 


118.457 


111.796 


108.522 


11 


171.358 


158.649 


136.252 


126.751 


117.556 


109.763 


105.366 


12 


171.589 


159.208 


136.927 


127.060 


117.121 


108.263 


102.728 


oo 


170.43(4) 


157.1(1) 


135.7(2) 


126.(3) 


115.1(1) 


103.7(1) 


87.6(7) 



Table 1. : Determination of the phase boundary between the high-temperature 
phase II and the IC-phase IV from the first crossing of the levels Oo and li. For 
brevity, we quote only the values for iV = 9, . . . , 12. However, in estimating the 
limiting values denoted "oo", data for all sites N = 3, . . . , 12 have been used. The 
estimated errors of the last given digit are quoted in brackets. 
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ip (degrees) 








N 


A =1.1 


1.2 


1.4 


1.8 


2.2 


2.6 


3.0 


9 


120.007 


124.409 


132.599 


143.899 


150.729 


155.270 


158 


526 


10 


116.541 


121.733 


131.092 


143.357 


150.491 


155.149 


158 


457 


11 


113.703 


119.695 


130.146 


143.192 


150.521 


155.232 


158 


553 


12 


111.355 


118.152 


129.613 


143.278 


150.718 


155.438 


158 


746 


oo 


99.(1) 


113.2(1) 


127.8(2) 


144.(1) 


151.(2) 


156.(2) 


159 


•(1) 



Table 2. : Determination of the phase boundary between the low-temperature 
phase I and the IC-phase III from the first crossing of the levels 1 and 2 2 . 
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Whether these transition lines have the properties of a Pokrovski-Talapov type 
transition |3^] is an interesting open question, because here two modulated phases 
have to coexist at the phase boundary line. 



3 Determination of the wave vector scaling expo- 
nent 

In this section we shall take a closer look at the iV-dependence of the level crossings 
in the low-temperature region A > 1. In Table 3 we give the cross-over positions in 
A for the superintegrable case which appear for iV up to 12 sites. 

This Table shows that the cross-over lines converge towards A = 1 and not towards 
the I/III-phase boundary A = 1.109518. So the low-lying p = 0-levels and their cross- 
overs do not care about the point where the IC-phase sets in. Between A = 0.901293 
and A = 1/0.901293 just the p = 0-levels are no longer the lowest levels. This is 
analogous to the formula ( |1.6| ) for the p = 0-gap valid for A < 1 which also takes no 
notice of the point A = 0.901293. 



N 


Q=0/l/2 1/0/2 1/2/0 2/1/0 2/0/1 


3 


1.8789156 








4 


1.3645610 








5 


1.2081036 


3.3430364 






6 


1.1368767 


2.0826986 






7 


1.0978185 


1.6742233 


4.8041526 




8 


1.0738838 


1.4761904 


2.8030552 




9 


1.0580721 


1.3610533 


2.1459063 


6.2689052 


10 


1.0470391 


1.2866459 


1.8229750 


3.5290973 


11 


1.0390113 


1.2350807 


1.6327694 


2.6244231 


12 


1.0329732 


1.1975179 


1.5083542 


2.1772063 4.255 


oo 


0.9997(5) 


0.999(1) 


0.990(6) 


0.93(5) 



Table 3. : Positions in A of the p = 0-sector lowest level crossings on the super- 
integrable line (p = 90° for periodic b.c. In the top line the order of the levels with 
different charges Q in the A-interval between the cross-overs is given. So the first 
column of A- values gives the position of the 0o/lo-cross-over, the second column the 
0o/2o-cross-over in the first excited state, etc. 

Critical exponents of the translationally invariant part of the spectrum at A = 1 
have been discussed in |29|, |24J. These authors find that scaling is asymmetrical 
around A = 1 and derived the specific heat exponent to be a = 1 — 2/n for the Z n - 
superintegrable models. From ( |1.6| ) it is immediate that the magnetic critical index 
is v — 1. The spacial correlation exponent v s has been calculated by Albertini et al. 
29| from the A- and iV-dependence of the Q = 0,p = 0-" ground" -state alone to be 
v s = 2/n. Baxter [23] finds from the exact solution of the superintegrable model the 
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same value 2/n for the interfacial tension exponent which desribes the exponential 
approach in N to degeneracy of the p = 0-lowest levels in the low-temperature phase 
near A = 1. 

We want to show that in the Z 3 -case ( |1.2| ) one obtains v s = | for the wave vector 
exponent at A = 1. We are not restricted by our method to the superintegrable point 
and find that v s varies little or not at all for a wide range of the chiral angle ip. All our 
results are consistent with the hypothesis that it is the inverse wave vector which sets 
the dominant spacial scale length of the p = sector of the spectrum around A = 1. 

The finite-size scaling rules for quantum chains with oscillating phases have been 



formulated in |33j . Let us define v s as the power of divergence of the inverse wave 
vector as A — > A c : 

l/K~(\-\ c )->". (3.1) 

In the region dominated by the length scale 1/K we use the scaling variable 

z = N 1 l v °{\-\ c ) (3.2) 

From Table 3 we deduce that A c = 1. In the scaling regime we can express the 
cross-over positions \ k (N) {k = 1, 2, . . .) in terms of values z k of the scaling variable. 
Calling Afc 5 fc/(iV) the distance between two level crossing points (e.g. at fixed ip) we 
can write 

A kyk ,(N) = A fc (JV) - \ k >(N) 

= {z k -z kl )N- 1 ^. (3.3) 

Using eq. ( |3.3| ) for iV and N — 1 sites, we get iVth approximants for v s : 

„ W »,fl = ™^» (3,, 

If our assumptions are consistent, the large iV-limit of these approximants has to give 
the same value v s for every choice of pciirs k^ k \ 

v s = lim u s (N,k,k'). (3.5) 

N— «x> 



One may also refer the distances to A c = 1 and use Afc iO (-^0 = Afc(iV) — A c in 
corresponding to z = in (|3~3|). 

In Table 4 we give numerical results from three different choices of Observe 
that from if = 120° down to <p = 60° we find no appreciable variation of v s with 
ip. In particular, there is no indication of a rise towards the ip = 0-Potts value 
v s = |. Unfortunately, for small (p the cross-over lines move to large A for our small 

available, so that the convergence of the u s (N, k, k') becomes poor e.g. for ip = 30°. 

From the level crossing points and the slopes of the gaps at these points, we 
can also obtain information on the specific heat exponent a and the magnetic gap 
exponent v [Q. The specific heat C v and its exponent a are defined by 

^(A) = ~^HA C -Ar. (3.6) 
In the region with ground-state cross-overs, its finite-size approximants C v (\, N) are 
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4S9^ 


4758 
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S840 


4S54 


5096 


^71 ^ 


9 


0.4961 


0.5050 


0.4891 


0.2636 


0.4255 


0.4719 


0.5489 


0.4154 


10 


0.5141 


0.5165 


0.5000 


0.3232 


0.4565 


0.5000 


0.5773 


0.4490 


11 


0.5276 


0.5252 


0.5093 


0.3682 


0.4806 


0.5223 


0.5977 


0.4753 


12 


0.5380 


0.5321 


0.5174 


0.4032 


0.4998 


0.5403 


0.6125 


0.4964 


oo 


0.59(1) 


0.59(2) 


0.66(4) 


0.61(2) 


0.67(4) 


0.66(2) 


0.66(1) 


0.66(1) 



Table 4. : Approximants v kjk '(N) to the wave- number exponent v s using three 
different levels choices of A k>k >: First A 10 , which is the distance of the first ground- 
state cross-over from A = 1, then A 2i o, the distance of the first cross-over in the first 
excited state from A = 1 and, in the last column, the distance A 21 between the first 
cross-overs in the ground state and first excited state. 



given by a sum of 5-function contributions (neglecting some smooth background): 

(3.7) 



fc-rnax (-^0 

C V (X,N) = Yl C k (N)6(X-X k (N)) 

k=l 



The weight factors of the 5-terms are given by the A-derivatives of the magnetic gaps 
at the ground state cross-over positions: 



C k (N) 



N 



d(E Qi - E Q . 



(3.8) 



A fc (A0 



The relevant charge sectors Q iy Qj rotate as A increases, see Table 3. 

In the scaling regime we express C V (X,N) in terms of the specific heat scaling 
function C(z) 

C v (X,N)=N a ^C(z) = N a ^Y,CkS(z-z k ), (3.9) 

so that 

C k (N) = C k N^~ l ^ Us . (3.10) 
From ( |3.10|) we get approximants a(N, k) for a: 



a(N, k) - 1 = u,(N, k) In 



'C k (N-l)\ (N-l 
/In 



C k (N) 



V N 



(3.11) 



Formula ( |3.11|) gives a practicable method for calculating a from the v s (N, k) and the 
(AQ 7^ 0)-magnetic gap slopes 
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The fact that the specific heat contributions from the ground-state crossings are 
proportional to magnetic gap slopes gives rise to the generalized hyperscaling relation 

a + v, + u = 2, (3.12) 

which is easily seen to be valid even for the finite- iV-approximants to a, v s and v fl33"f . 

For detailed tables of the relevant sequences of approximants, we refer the reader 
to Here we just mention that we find v = 1 within a few percent error for all 

60° < ip < 120°. So for these values of ip eq.(glg) reads § + 1 + | = 2. For the Potts 
case ip = 0° we know that it must read | + | + | = 2, but we are not able to tell at 
which ip the first equation goes over to the second one. 



4 Particle interpretation of the low-lying spectrum 

In this section we attempt to describe the low-lying part of the spectrum of TC of 
eq.( |1.2[ ) in terms of elementary excitations. For the superintegrable case such a dis- 
cussion has been given e.g. in [[38, 26], but it is instructive to consider the dependence 
on the chiral angle too. We shall treat the self-dual case ip = <f> if not stated otherwise. 

Let us start with the information we have for ip = 0. This is the standard 
Z3-Potts model which has a second order phase transition at A = 1 decribed by a 
c = |-conformal field theory |39|, |D|. Zamolodchikov has shown |^1J that the high- 



temperature scaling regime around A = 1 is described by the thermal perturbation of 
the c = |- conformal theory, which, as he showed, is still integrable. The S"-matrix 
of this massive theory factorizes and contains just a single pair of particles with Z 3 - 
charges Q = l and Q = 2. Because of the additional ^-symmetry present at ip = the 
Q — 1 and Q = 2-sectors are degenerate and consequently both particles have the same 
mass, forming a particle-antiparticle pair. The higher levels of the high-temperature 
scaling Potts-model are multiparticle scattering states of which even the phase shifts 
are known. 

In order to find out whether the spectrum for ip ^ is continuously obtained 
from the Potts spectrum, we first take a look at a high-temperature expansion of the 
p = 0-states. In the interval — | < ip < n we find for the lowest Q = 1, p = 0-gap up 
to order A 2 : 

. _ — . / , N . . it - ip 4A ip 2A 2 sin|(7r — 2(p) „. . /A . 

AE 1>0 &, A) = 4 sin - -= cos \ - -= . \) / / (cp ) + . . . 4.1 

3 V3 3 V3 sm |(7r + 2ip) 

Here /(<£>) is a smoothly varying function of ip, which satisfies 

/(-f) = 2/ v / 3; /(0) = 1; /(f ) = V3/2; lim /(^)/sin|(7r - ip) = 2y/3. (4.2) 

Applying a CP-transformation, the same formula can be used for the Q = 2, p = 0-gap 
if \ip\ < f: 

AE 2fi (<p,\) = AE lt0 (-<p,\). (4.3) 
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For small values of A, for which these formulae are useful, the level Eo(Q = l,p = 0) 
is isolated in the spectrum of the Q = 1-charge sector if < p < tt. Similarly,the level 
E (Q = 2,p = 0) is isolated in the spectrum of the charge sector Q = 2 for < <p < |. 
So we tentatively call these two levels "single particle levels" and assign the masses 

mi = AE lfi ; m 2 = AE 2fi . (4.4) 

For <p — > both particles become degenerate, in agreement with our discussion of the 
Potts case. 

Now let us take a look at formulae (|4.1|) ( |4.3|) for <p — > |. Disregarding first the 
A 2 -term, with <p increasing, m 2 increases faster than mi decreases. For ip = | (the 
superintegrable case) we seem to get 

m l = 2-2A + 0(A 2 ); m 2 = 4 - 2A + 0(A 2 ). (4.5) 

While this formula is allright for m\ (we even know that all higher order terms in A 
vanish as p — > |-) it is not correct for m 2 , since for m 2 the coefficient of A 2 diverges 
as p — > |. Physically this divergence is not surprising because as p — > | particle m 2 
moves towards the threshold of the two-mi-particle scattering states which are in the 
same Q = 2-sector, since the Z3 charges of the particles are additive mod 3. 

Calculating numerically ni\ and m 2 for many points of the region < p < |, < 
A < 1 we find that there is a smooth change of the ratio m 2 /mi from its Potts- value 
m 2 /mi = 1 to m2 /mi — ► 2 as p approaches tt/2 from below. 

Since we can easily calculate numerically the about eight lowest p = levels of 
each charge sector for all < p < 7T, we now check whether the pattern of higher levels 
fits to the proposed particle interpretation. In Figs. 3 and 4 we show the expected 
patterns of levels and thresholds, which should be qualitatively different for p < | 
and p > tt if in the latter case m 2 becomes larger than 2mi. 

• Sector Q = 0: As long as m 2 < 2m x (p < |), the first threshold is due to 
the scattering of two different particles mi and m 2 , followed at AE = 3mi by 
the 3-mx-particle scattering states (here again we refer the energy gaps to the 
Q = 0, p = 0-ground state). This two-particle threshold moves into the 3mi- 
continuum at the superintegrable line if there we have m 2 = 2mi. In Table 5 we 
give the finite- N energy gaps accessible to our numerical diagonalization of 7H 
for the superintegrable case and A = 0.5. We choose this value of A because then 
from ( |1.6|) we have m x = 1 and a normalization of the data is unnecessary. The 
superintegrable case is particularly interesting, because it makes it necessary to 
disentangle physically different but degenerate levels. 

Not surprisingly after this general discussion, in our numerical calculation, the 
lowest bunch of levels appears at AE = 3mi. Now, in order to find out, which of 
these levels are mi + m 2 -states and which are the mi + mi + mi-states, we look 
into their finite-size behaviour. For a single particle level we expect exponential 
convergence for iV — ► 00: 

AE(N) - AE(oo) = exp(-JV/0 + • • • (4-6) 
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where £ is a correlation length. In contrast, for multiparticle scattering states 
we expect power convergence |42|, [43| : 



AE(N) - AE(oo) ~ N~ y + ... (4.7) 
We calculate approximants y n to y using neighbouring values of N: 

/ AE(jV) - AE(oo) \ / JV_v 

yjv V As(iv - 1) - as(oo) y 7 ViV-iy v ; 

Table 6 shows some y ra calculated from the iV-dependence of the energy levels 
in Tab. 5. We see that the limiting values of y come out clearly to be very 
close either to y = 2 or y = 3 and allow us to decide which of the levels are 
mi + m 2 -states and which are m\ + m x + m\. 

• Sector Q — 1: Here we see exponential convergence for the isolated level mi, and 
with good precision we see at AE = 4 the bunch of the next higher levels. There 
are levels converging with y = 2 and others converging with y = 3. The obvious 
interpretation is in terms of m,2 + va-i- two-particle states and of ni\ + mi + mr 
three particle states. We have checked that moving towards smaller values of <p 
these levels split in the expected way: for m 2 < 1m\ the levels interpreted as 
"two-particle" go lower than the assumed three-particle states. 

• Sector Q = 2: For e.g. if = 45° there is the clear exponential convergence of 
the isolated m2-level. However, the interesting and crucial test comes now in 
the superintegrable case: Among the several levels at AE = 2 (we still consider 
A = 0.5) which show power convergence, can we distinguish one level which 
shows exponential behaviour? Already without calculating the sequences y n 
we see that the level in the first row of the bottom part of Table 5 sticks out 
by its very fast convergence. This is substantiated in Table 6 where we see 
the corresponding sequence of the y n growing very fast. So we have evidence 
that the particle m 2 survives up to ip = |. We have tried to follow into 
the continuum for (p > |. For (p = 105° and A = 0.15 we have succeeded in 
distingishing one level which converges faster than a low power y = 2 or y = 3. 

We have also performed calculations for <p = 45°, 105°, both for the self-dual 
and for the integrable model. In all these cases we find a pattern of the ca.6 
lowest levels of each charge sector as scetched in Figs. 3 and 4. 

By analogy, we expect that the superintegrable Z n -models have low-lying excita- 
tions corresponding to n — 1 particle species, one in each Q ^ 0-sector with mass 
m Q = Qm x . 

Fig. 3: (following page) General structure of the low- lying spectrum. Left hand 
side of the picture: Structure of the spectrum in the different charge sectors if the chiral 
angle <p has a value below (p — |. Right hand side: for p above the superintegrable 
value. 
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Table 5: The lowest energy gaps AEqj, as defined in eq.( |1.5|) , for the Z 3 -hamil- 
tonian eq.( [1.2|) , for the superintegrable case (f> = <p — ir/2, and A = 0.50. Upper 
Table: sector Q = 0, middle table: sector Q = 1 and bottom table: sector Q = 2. 
The numbers given in brackets indicate the estimated error in the last written digit. 
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Table 6: Exponents i/n of the convergence iV — > oo as defined in eq. (|4.8|) , for the 
superintegrable case and A = 0.50 for the lower levels given in Table 5. "expon." in 
the bottom row means that the increasing sequence of the in the column above 
indicates exponential convergence. 



Having been successful with the two-particle interpretation of the low-lying spec- 
trum, we finally consider the energy-momentum relation of these particles. Since 
parity is generally not a good quantum number of our model, the dispersion relation 



is not expected to be symmetrical between positive and negative momenta. In (19 
such unsymmetrical dispersion curves have been given for the superintegrable case. 
We have studied the (^-dependence of this asymmetry. Of course, for (p —>■ the 
asymmetry must vanish because for (p = <ft = the hamiltonian 7i is parity invariant. 
The high-temperature expansion gives us some information about the increase of the 
asymmetry with ip (we consider again the self-dual case): 

For calculating the limit N — > oo using lattice momenta p as in eq.( |1.4| ), we 
introduce the macroscopic momentum P defined by 

P =-^P (-vt>P>tt) (4.9) 



Then, in generalisation of (4J.), (43) we find up to order A 



AE 1>0 (\, P) = 4 sin - i| cos(P -\<p) + ... (0 < <p < tt) (4.10) 

and 

AP 2i0 (A,P) = 4sin^^---^cos(P + ^) + ... (0 < <p < tt/2). (4.11) 

3 v 3 

This shows that for small A, the asymmetry of the energy-momentum relation with 
respect to P — ► — P is due only to the presence of a macroscopic momentum P m = 
±y?/3 (observe that this is the same momentum which appears in eq. (p.l|) for (p = 
|). We have checked numerically that for A < 0.1 these formulae are an excellent 
description of the data. Only for large values of A the situation is getting complicated. 
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Of course, since we considered just a few low-lying levels, we are not able to make 
statements about higher levels. In ref . 1 26| it has been shown that in the superinte- 
grable case the model contains levels which cannot be interpreted as quasiparticles in 
the sense that their energies and momenta receive unrelated additive contributions. 



5 Conclusions 

In this talk we reported results of a numerical investigation of the low-lying spectrum 
of the Z3-chiral Potts hamiltonian ( |1.2|) for N = 2, . . . , 12 sites. We get an overview 
about the dependence of ca. 18 levels on the chiral angle and the inverse temperature. 
This fits to the interpretation of the high-temperature low-lying spectrum in terms of 
two basic particles with Z 3 -charges Q — 1 and Q = 2 and their scattering states. The 
superintegrable case is distinguished by the special 1 : 2-mass ratio of the particles. 
Even when moving slightly into the continuum formed by two Q = 1-particles, the 
heavy Q = 2-particle can be distinguished by its exponential finite-size behaviour. 
Two- and three-particle scattering states can be distinguished by their different finite- 
size power behaviour. High-temperature perturbation theory gives a simple expression 
for the effect of parity violation on the energy-momentum relation of the particles. 

Particular attention has been devoted to the study of ground-state level crossings 
in the low-temperature regime. We found that the massive phases have oscillatory 
correlation functions and calculated the wave-vector critical exponent from the low- 
temperature phase. It turns out that for translation invariant features of the model 
in the neighbourhood of the self-dual line, the diverging inverse wave vector defines 
the most relevant length scale. 
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